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Abstract 

We study the influence of a possible coherent component in the 
boson source on the two-, three- and n-particle correlation func- 
tions in a generalized core/halo type of boson-emitting source. In 
particular, a simple formula is presented for the strengh of the 
n-particle correlation functions for such systems. Graph rules are 
obtained to evaluate the correlation functions of arbitrary high 
order. The importance of experimental determination of the 4-th 
and 5-th order Bose-Einstein correlation function is emphasized. 



1 Introduction 



Intensity correlations were discovered first in astrophysics by R. Hanbury Brown 
and R. Q. Twiss who invented this method to determine the angular diameter 
of main sequence stars (HBT effect). In particle physics, intensity correlations of 
pions were observed by Goldhaber, Goldhaber, Lee and Pais (GGLP effect) 
Bose-Einstein correlations are intensity correlations among detected bosons, that 
are studied mainly with the purpose of reconstructing the space-time picture of 
particle production. The analysis of higher-order Bose-Einstein correlation func- 
tions became a focal point of current research interest. 

In particle physics, significant three or higher order Bose-Einstein correlations 
have been extracted from the data sampled by the AFS f|, the NA22 (|, §, [| 
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and the UAl-collaborations M. These data were used to test the possible exis- 
tence of a coherent source in multi-particle physics and to compare the correlation 
functions to the strength of these correlations predicted from the quantum optical 
QO formalism g g, As the precision of the measurements improved, the QO 
predictions with higher order correlations were found to be marginally consistent 
with the data on 3 and 4 order Bose-Einstein correlation functions (BECF-s) 
in (tt + /K + ) +p reactions at CERN SPS. Recently, this basic QO formalism was 
shown be insufficient to simultaneously describe the high precision UA1 data on 
two- and three-particle Bose-Einstein correlations [T^| . 

In high energy heavy ion physics, the first experimental determination of the 
three-particle correlation function has just been reported by the NA44 collabora- 
tion H, 0. NA44 reports that the genuine three-particle correlation is 
quite suppressed in the studied reaction, S + Pb collisions. By the genuine three- 
particle correlation is meant the part of the three-particle correlation that is not 
due to included combinations of two-particle correlations. This suppression can 
be expressed as a phase factor, cos(0), of the three-particle correlation function in 
the case of totally incoherent particle production. In that case this phase factor is 
related to an asymmetry of the particle source, not possible to extract from two- 
particle correlations. Theoretical estimates of this asymmetry effect on the phase 
factor show very small departures from cos(^>) « 1. [|l7| [l8|, [TjJ. The large departure 
from cos(</>) = 1 found by the NA44 collaboration, cos(</>) = 0.2 ± 0.2 pi, ought 
to be due to some other mechanism. We will discuss the possibility of a partial 
coherent source in this Letter. A possible existence of such an extra phase in the 
three and higher order correlation functions was noted already e.g. in papers by the 
NA22 collaboration but no experimental evidence has been put forward for a 
cos(</>) ^ 1 value in particle physics. 

From the theoretical side, Cramer and Kadija predicted up to order 6 the 
strength of Bose-Einstein correlations for sources with partially coherent and inco- 
herent components that included also a possible contamination by mis-indentified, 
non- interfering particles [ pp| . Their formulas were obtained in the quantum-optical 
formalism. Recently, Suzuki and collaborators calculated higher order exclusive 
Bose-Einstein correlations from the generating functional approach to quantum- 
optical formalism for the case that the source has M incoherent and one coher- 
ent component. 

Recently, multi-particle symmetrizations up to arbitrary high order were evalu- 
ated exactly by Zhang |2^| for the special case of a pion-laser model proposed by 
Pratt in ref. p3f . Surprizingly, the structure of the n-particle inclusive correlation 
functions in terms of the Fourier-transformed inclusive emission function was found 
to be the same as the structure of the n-particle exclusive correlation functions in 
terms of the single- particle exclusive emission function |p2[ . However, this result 
is valid only in case when Bose-Einstein condensation, hence the development of 
partial coherence, is not yet reached p4| . 

A simple recurrence relation was obtained for the strength of the higher order 
correlation functions of core/halo type systems p5|. Such systems are boson emit- 
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ting sources where some particles come from the a incoherent center of the particle 
emission, that is assumed to be resolvable by the Bose-Einstein microscope. The 
rest of the particles is assumed to come from the halo region, that corresponds 
to large length-scales not resolvable by intensity interferometry |2(| ^7|. In Let- 
ter |2jj] , a prediction was made for the strength of third order and arbitrary order 
BECF assuming that the core has no coherent component. 

The purpose of the present Letter is to investigate the effect of a partially co- 
herent component in the core of particle emission. We present a generalization of 
the earlier recurrence relations in ref. p5|, the new expressions also yield an easy 
way to calculate formula for the strength of the n-th order correlation function with 
a partially coherent and a halo component and we apply these expressions to the 
NA44 data on S + Pb collisions. 



2 Basic definitions 

The central assumption of the core/halo model is that the reduction of the intercept 
parameter of the n-particle BECF-s is due only to the presence of the long-lived 
resonances j25|. This assumption was motivated by the success of fully incoherent 
event generators like RQMD or VENUS in the description of two-particle BECF-s. 

The emission function of the whole source can be written as a sum of a contri- 
bution from the core and from the halo, where halo stands for the decay products 
of the (unresolvable) long-lived resonances. The core is indexed with (c) , the halo 
is by (h). 

S(x, k) = S c (x, k) + Sh(x, k) (1) 

In earlier studies of the core/halo model it was assumed that S c (x, k) describes 
a fully incoherent (thermal) source. Now we assume, that some fraction of the core 
emits bosons in a coherent manner, e.g., due to emerging formation of pion lasers or 
Bose-Einstein condensates of pions or production of disoriented chiral condensates 
or so we define 



S c (x,k) = S*(x,k) + Si(x,k) (2) 

where the upper index p stands for coherent component (p as partial), upper index 
i stands for incoherent component of the source. 
The invariant spectrum is given by 

N(k) = J d 4 xS(x,k) = N c (k)+N h (k) (3) 

and the core contribution is a sum of the coherent and incoherent components: 

N c (k) = f d 4 xS c (x, k) = NP(k) + N l c (k) (4) 
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One can introduce the momentum dependent core fractions / c (k) and partially 
coherent core fractions p c (k) as 

/ c (k) = N c (k)/N(k) (5) 
p c (k) = iV c p (k)/iV c (k) (6) 

The halo and the incoherent fractions fh, fi are 

/ & (k) = N h (k)/N(k) = 1 - / c (k) (7) 
/,(k) = iVj(k)/7V c (k) = l-p c (k) (8) 

Note that our definition of the momentum dependent, partially coherent core 
fraction, p c (\t) should be clearly distinguished from the chaociticy p of Weiner f2§| , 
defined as p — (n c hao) I ("-tot) , the ratio of the mean number of particles from the 
chaotic source to the mean total multiplicity. If we neglect the momentum depen- 
dence of / c (k) and p c (k), the core fraction and the partially coherent core frac- 
tion, formally one obtains p = 1 — p c f c . However, we distinquish the resolvable 
intercept A* from the exact intercept A KCt , in contrast to ref. pq| . For exam- 
ple, in case of two-particle correlations, A*,2 = /<?[(! — Pc) 2 + 2p c (l — Pc)], while 
^xct,2 = A».2 + (1 — f c ) 2 + 2/ c (l — f c ) hi our case, while in case of the quantum 
optical formalism without long lived resonances, A^° = 2p(l —p) +p 2 = 1 — (1 —p) 2 . 



3 The strength of the n-particle correlations, A* n 

We define the n-particle correlation function as 

C„(l,2,.-,n) = C n (k 1 ,k 2 ,...,k»)= (9) 

7Vi(ki)iVi(k2)...7Vi(k„) 

iV n (l,2,...,n) , s 

(10) 



N 1 (l)N 1 (2)...N 1 (n) 



where a symbolic notation for k^ is introduced, only the index of k is written out in 
the argument. In the forthcoming, we shall apply this notation consistently for the 
arguments of various functions of the momenta, i.e., /(ki, k, , k m ) is symbolically 
denoted by f(i, j, m). 

We find that the intercept of the n-particle correlation function ( extrapolated 
from finite relative momenta to zero relative momentum) is given by the following 
formula, 



c n (h = k h Vij) = i + A., n = 1 + E ( f) a jfi [(! -Pc) j +m(i -PcY- 1 ] 



(11) 
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where ay counts the number of fully mixing permutations of j elements. This can 
be calculated from a simple recurrence, as obtained in ref . Q . 

Note that the equations of ref. |2(| were given for the purely incoherent core, 
and they are modified above for an additional coherent component in a straight- 
forward manner. In general, terms proportional to f£ in the incoherent case shall 
pick up an additional factor [(1 — p c ) J + jp c (l — PcY^ 1 ] m case the core has a 
coherent component. This extra factor means that either all j particles must come 
from the incoherent part of the core, or one of them can come from the coherent, 
the remaining j — 1 particles from the incoherent part. If two or more particles come 
from the coherent component of the core, the contribution to intensity correlations 
vanishes as the intensity correlator for two coherent particles is zero. 

Let us indicate the number of permutations that completely mix exactly j non- 
identical elements by ay . There are exactly f different ways to choose j different 

elements from among n different elements. Since all the n\ permutations can be 
written as a sum over the fully mixing permutations, the counting rule yields a 
recurrence relation for Oj, ref. |2q]: 

«n = nl-l-X;^)^. (12) 

The first few values of otj are given as 

a t = 0, (13) 

a 2 = 1, (14) 

as = 2, (15) 

a 4 = 9, (16) 

a 5 = 44, (17) 

a 6 = 265. (18) 

We have the following explicit expressions for the first few intercept parameters: 

A*, 2 = / c 2 [(l- Pc ) 2 + 2p c (l- Pc )] (19) 

A*, 3 = 3/ c 2 [(l-p c ) 2 + 2p c (l- Pc )] 

+2f'![(l-Pcf + 3 Pc (l- Pc ) 2 } (20) 

A*, 4 = 6/ 2 [(l-p c ) 2 + 2p c (l-p c )] 

+8f?[(l- Pc ) 3 + 3 Pc (l- Pc ) 2 } 
+9/ c 4 [(l-p c ) 4 +4 Pc (l- Pc ) 3 ] (21) 

A*, 5 = 10/ 2 [(l- Pc ) 2 + 2p c (l-p c )] 

+20/' 3 [(l-p c ) 3 +3 Pc (l-p c ) 2 ] 
+45/ c 4 [(l-p c ) 4 + 4 Pc (l-p c ) 3 ] 
+44/ c 5 [(l-p c ) 5 + 5p c (l-p c ) 4 ] (22) 
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In the above equations, the effective intercept parameters, the core fraction and the 
partially coherent fraction are evaluated at a mean momentum K, A*,„ = A* in (K), 
f c = / C (K) and p c = p c (K). 



4 The n-body correlation function 

Let us give the closed form for the full correlation function for arbitrary high order of 
correlation functions, generalizing the results of ref. |^5) for an additional partially 
coherent component in the source: 

Let stand for those permutations of (1, ri) that are mixing all the numbers 
from 1 to n and let us indicate by pi the value which is replaced by i in a given 
permutation belonging to the set of permutations pW. (Superscript indexes a set 
of permutations, subscript stands for a given value). Then we have Pi ^ i for all 
values of i — 1, n. 

If the partial coherent component is vanishing, the general expression for the 
n-particle inclusive correlation function C n (l, n) was given in ref |2q | as 

n n / j 

C n (l,...,n) = 1 + E E En^^'VJ- (23) 

j=2 ix,...,ij=l p(j) k=l 

Here ^ indicates that the summation should be taken over those set of values of the 
indices which do not contain any value more than once, and the Fourier-transformed 
emission function of the core is 

ic(i,i) = ||| (24) 



h(iJ) = KU,i)-gfr^*r B U,i), (25) 
b c (i,i) 

In the above equations, the tilde denotes Fourier-transformation over the relative 
momenta, 

S c {l,m) — J d 4 xexp[i(ki - k m ) ■ x]S c (x, — + ^ m ) (26) 

and similar expressions hold for the coherent and the incoherent components of the 
core. [] 

The expression in eq. ( p3| ) is valid not only for the case when exactly n bosons 
are in the system, full symmetrization is performed, C n (l, 2, n) stands for the n- 
particle exclusive correlation function and S c (i,j) stands for the Fourier-transformed 

1 Note that with this definition the normalized Fourier-transformed emission function becomes 
asymmetric to the exchange of the arguments and complex conjugation: although the relationship 
S c (i,j) = is satisfied. 
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core emission function without modifications due multi-particle symmetrization j22|, 
p4j . In addition, eq. (^3|) is also valid when the only source of correlations between 
pions is due to Bose-Einstein symmetrization, the number of pions is randomly vary- 
ing from event to event, and C n (l, 2, n) is interpreted as the n-particle inclusive 
correlation function p4 | , and S c (i,j) includes all higher order symmetrization 
effects. However, eq. ([23D is valid only if the core has no partially coherent com- 
ponent. If a coherent component is present, one can introduce the normalized 
incoherent and partially coherent core fractions as 

= |m (27) 

and we obtain 

n n i ( j 

C„(l,...,?i) = 1 + ^2 H \ II M m k)[l -Pc(m k )} s l c (m k ,m Pk ) 

j=2 mi,...,mj=l nU) Kk=l 
j j 1 

+ ^2 fc(mi)Pc(mi) s p c {m h m Pl ) JJ f c (m k )[l - p c (m k )] s l c {m k , m Pk ) \ (29) 
(=1 fe=l,fe^i J 

This expression contains phases in the Fourier-transformed, normalized source dis- 
tributions. Actually, two (momentum dependent) phases are present: one denoted 
by z (k m , k„) in the Fourier-transformed normalized incoherent core emission func- 
tion, s l c (k m , k„) and another independent phase denoted by </> c (k m , k„) is present in 
the the Fourier-transformed normalized coherent core emission function, s£(k m , k n ). 
One can write 

S 4 c (k m ,k„) = |S*(k m ,k„)|exp[^(k m ,k n )], (30) 
S p c (k m ,k n ) = |SP(k m ,k n )|exp[^(k m ,k„)]. (31) 

The shape of both the coherent and the incoherent components is arbitrary in these 
equations, but should correspond to the space-time distribution of particle pro- 
duction. If the variances of the core are finite, the emission functions are usually 
parameterized by Gaussians. If the core distributions have power-law like tails, like 



in case of the Lorentzian distribution |29|, then the Fourier-transformed emission 



functions correspond to exponentials or to power-law structures 1 30 . For complete- 
ness, we list these possibilities below: 

|S*(k m ,k„)| 2 = eM-R*Qln) or (32) 
\s l c (k m ,k n )\ 2 = exp(-RiQ m n) or (33) 
|^(k m ,k„)| 2 = ^(iWmn) 6 * etc., (34) 
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fc 


P 


A*,2 


A*, 3 


A M 


A*, 5 


0.60 


0.00 


0.36 


1.5 


5.1 


17.2 


0.70 


0.50 


0.37 


1.4 


4.3 


11.9 


1.00 


0.75 


0.44 


1.6 


4.3 


10.5 



Table 1: Evaluation of the strength of higher order correlation functions, A*, ra , for 
various core fractions and partially coherent fractions allowed by NA44 two- and 
three-particle correlation data. 



|S?(k m ,k„)| 2 = exp(-i^Q^„) or (35) 
|s£(k m ,k„)| 2 = exp(-R p Q mn ) or (36) 
|^(k m ,k„)| 2 = a p (R p Q mn p etc.... (37) 

In the above equations, subscripts i and p index the parameters belonging to the 
incoherent or to the partially coherent components of the core, and Q mn stands for 
certain experimentally defined relative momentum component determined from k m 
and k„. 

A straightforward counting yields that in the limiting case when all momenta 
are equal, the simple formula of eq. ( |ll| ) follows from the shape of the n-particle 
Bose-Einstein correlation functions of eq. as s c (i,i) — s^(i,i) = 1. 



5 Application to three-particle correlation data 

As an application of the above formalism, we attempt to determine the core frac- 
tion f c and the partially coherent fraction p from the strength of the NA44 two 
- and three- particle correlation functions, A,^ and A*^, in the CERN SPS S + 
Pb reactions. The two experimentally determined values are A*,2 = 0.44 ± 0.04 
and A».3 = 1.35 ± 0.12 (statistical errors only). Figure 1 illustrates the 2 a 
contour plots in the (f c ,Pc) plane, that is obtained for these parameters from the 
experimental values of A*^ and A*^. 

The overlap area in Fig. 1 shows, that a big range of (f c ,Pc) values is able 
to describe simultaneously the strength of the two-particle and the three-particle 
correlation functions within two standard deviations of experimental errors. Thus 
neither the fully incoherent, nor the partially coherent source picture can be ex- 
cluded. 

Now we can predict the intercept of higher order correlations to see if they 
become more sensitive to the presence a partial coherent source. In Table 1 we have 

2 Coulomb corrections are large in heavy ion collisions and the value of A*, 3 was determined with 
the help of a newly developed Coulomb 3-particle wave-function integration method described in 
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0.2 0.4 0.6 0.8 1 



Pc 

Figure 1: Allowed regions for possible values of the f c core fraction and the p c 
partially coherent fraction are evaluated on the two a level from the intercept of 
the second order and the third order correlation functions, A* ; 2 and A,^. 
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evaluated the the A*, 2, A*,3, A*, 4, A^s-values for some cases in the overlap region. 
We find that the A*,5 is almost a factor of 2 larger for a completely incoherent 
source, than for a partially coherent source with no halo component, although 
within experimental errors both cases describe A* : 2 and A*,3. This is in agreement 
with Cramer and Kadija, who have pointed out that for higher values of n the 
difference between a partially coherent source and the fully incoherent source will 
become larger and larger Q ■ 

The results presented here imply that the measurement of higher order cor- 
relations, 5-th order, is necessary to determine the value of the degree of partial 
coherence of the source in this reaction. 

6 Summary, conclusions 

In summary, we have found a simple generalization of the core-halo model for the 
case when the core has a partially coherent component. The strength of the n- 
particle correlation function can be evaluated for arbitrary value of n with the help 
of a simple recurrence formula. 

The shape of the n-particle Bose-Einstein correlation function was determined 
in terms of the Fourier-transformed emission function of the incoherent and the 
partially coherent component of the source. The graph rules for the calculation of 
these functions are summarized and illustrated graphically in Appendix A. 

We found that the strengths of the second and the third order Bose-Einstein 
correlation functions in the NA44 S + Pb reaction at CERN SPS can be accommo- 
dated simultaneously both in a fully incoherent core picture (p c = 0) with a halo 
fraction of f c — 0.6 as well as in a partially coherent core picture that has no halo 
component, p c — 0.75 and f c = 1. However, the strength of the fourth and fifth 
order correlation functions is shown to be quite different in the two scenarios. 
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Appendix: Graph rules 



A straight-forward calculation of the higher order Bose-Einstein correlations for a 
partially coherent core/halo type of systems is possible with the help of a set of 
graph rules that we determine below. Although the graphs we describe are similar 
to those of ref. J2l| , the rules are different as we have multiplicative factors for each 
vertex (that carry one momentum label each) and for each line (that connect two 
vertex, hence carry two momentum labels each) . 

Figures 2 and 3 graphically illustrate the rules of calculations of the contribu- 
tions of the incoherent and coherent core components to the n-particle correlation 
function C„(ki, k„) for the cases n = 2, 3 and 4, respectively. Circles can be 
either open or full. Each circle carries one label (e.g. j) standing for a particle with 
momentum kj . Full circles represent the incoherent core component by a factor 
/c(j)[l ~ Pc{j)]i whereas open circles correspond to the coherent component of the 
core, a factor of f c {j)p c (j), as defined in eq. <Js| — [sj) and also shown on Fig. 2. 

For the n-particle correlation function, all possible j-tuples of particles have to 



be found. Such j-tuples can be chosen in . I different manner. In such a j-tuple, 



either each circle is filled, or the circle with index k is open and the other j — 1 circle 
is filled, which gives j + 1 different possibilities. All of the permutations that fully 
mix either j = 2, or 3, ... , or n different elements have to be taken into account for 
each choice of filling the circles. The number of different fully mixing permutations 
that permute the elements i±, ...ij is given by ccj and can be determined from the 
recurrence of eq. (|l^). 

Lines connect two circles (or vertexes) denoted e.g. by The lines stand 

for factors that depend both on k^ and on kj. Full lines represent incoherent 
- incoherent particle pairs, and corresponds to a factor of s l c (i,j). Dashed lines 
correspond to incoherent-coherent pairs, and carry a factor of s^(i,j). The lines are 
oriented, they point from circle i to circle j, corresponding to the given permutation, 
that replaces element j by element i. Dashed lines must start from an open circle 
and point to a full circle. 

All possible graphs must be drawn that are in agreement with the above rules. 
The result corresponds to the fully mixing permutations of all possible j-tuples 
(J = 2, ...n) chosen in all possible manner from elements (1, 2, n). 

Each graph adds one term to the correlation function, given by the product of all 
the factors represented by the cirles and lines of the graph. Note that the directions 
of the arrows matter, as reflected by the unequality in eq. (|2^). The correlation 
function C(l, ...,n) is given by 1 plus the sum of all the graphs. 

Finally we note, that for the n-particle cumulant correlation function, n circles, 
representing the n particles, should be connected in all possible manner correspond- 
ing only to the fully mixing permutations of elements (1, n). Disconnected graphs 
do not contribute to the cumulant correlation functions, as they correspond to per- 
mutations, that either do not mix all of the n elements or can be built up from two 
or more independent permutations of certain sub-samples of elements (1,2, n). 
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Figure 2: Graphs determining the second and the third order correlation function 
for partially coherent core/halo sources. 
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